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Abstract 



We study classical solutions of a particular version of the modified Skyrme 
model in (3+1) dimensions. The model possesses Skyrmion solutions as well as 
stable domain walls that connect different vacua of the theory. We show that 
there is an attractive interaction between Skyrmions and domain walls. Thus 
Skyrmions can be captured by the domain walls. We show also that, when the 
mass term is of a special type, the model possesses bound states of Skyrmions 
and of the domain wall. They look like deformed 2-dimensional Skyrmions 
captured by the wall. The field configurations of these solutions can inter- 
preted as having come from the evolution of the 3-dimensional Skyrmions 
captured by the domain wall. For more conventional choices of the mass 
term of the model in the model the attraction between the Skyrmions and 
the wall leads to the capture of the Skyrmions which are then turned into 
topological waves which spread out on the wall. We have observed, numer- 
ically, such captures and the emission of the waves. We speculate that this 
observation may be useful in the explanation of the problem of baryogenesis 
and baryon-antibaryon asymmetry of the Universe. 
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I. INTRODUCTION 



The Skyrme model in (3+1) dimensions, originally introduced by Skyrme [I]], has at- 
tracted a lot of attention as a reasonable good effective theory for describing baryons 
@. Baryons in this theory are described by topologically nontrivial classical configura- 
tions (Skyrmions). Recently there has been some discussion of the problem of the baryon- 
antibaryon asymmetry of the Universe. This problem may be related to the problem of 
interaction between Skyrmions (or baryons) and domain walls during the early stages of the 
Universe. 

The interaction of particles with domain walls has become a subject of interest following 
a seminal paper by Voloshin ||. A detailed analysis of the scattering of Abelian gauge 
particles on domain walls may be found in a paper by Farrar et al ||]. The problem of the 
chiral fermion determinant in the presence of a domain wall was discussed in ||. It is also 
worth mentioning that the expanding bubbles of a new phase during the electroweak phase 
transition may be an additional source for the CP - violationg effects of the electroweak 
baryogenesis, see, e.g. f|0]. 

Recent interest in the problem of the interaction of particles with domain walls was 
generated by a paper of Dvali et. al [Q. In this paper the authors consider the interaction 
of monopoles with domain walls. They suggest that this interaction may help to solve the 
problem of low density of monopoles in the Universe. They conjecture that the domain walls 
would sweep up all the monopoles which would then defuse along the domain walls. 

Very recently this conjectured interaction of monopoles with domain walls was supported 
by numerical stimulations of monopole-domain wall collision processes as reported in |J. 

The fact that extended objects, such as Skyrmions, interact nontrivially with domain 



walls had already been discussed in our previous paper JT0[. In that paper we studied the 



interaction of Skyrmions with domain walls in a two-dimensional baby-Skyrmion model |TT 



In |10| we showed that domain walls can attract and absorb Skyrmions. After its collision 
with a domain wall the Skyrmion splits into two parts which propagate (with the speed of 
light) along the domain wall. Each part of the former Skyrmion looks like a wave packet 
which carries one half of the topological charge of the Skyrmion. That is why the compound 
objects like Skyrmions and monopoles can interact nontrivially with domain walls. This 
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interaction requires further careful study. 

In this paper we concentrate our attention on the problem of the interaction of Skyrmions 
with domain walls in a realistic 3-dimensional case. To do this we consider the conventional 
Skyrme model |],[| but with a slightly modified expression for the mass term. Our modifi- 
cation is motivated by the existence of solutions in the form of domain walls. In fact, first 
we introduce the simplest modification which is then further altered by the introduction 
of an additional parameter. The model is introduced and discussed in detail in the next 
section. For our choice of the mass term it possesses solutions which correspond to both 
Skyrmions and domain walls. Next we study the problem of their interaction. We demon- 
strate that there exists a solution, which may be considered as a bound state of a Skyrmion 
and a domain wall. We also show that the Skyrmions and the domain walls attract each 
other. We have performed a series of numerical simulations studying this interaction. In our 
simulations we have always found the attraction although the final stage of the evolution 
process(ze whether the topological waves on the wall spread out or oscillate) depends on the 
details of the modified mass term of the model. In each case the Skyrmions are captured, 
resulting in the emission of topological waves. The waves spread out on the wall in the 
model with the simplest version of the modified mass term, or oscillate when the mass term 
is further modified to prevent this spreading. 

II. THE MODEL 

To study the interaction between Skyrmions and domain walls in (3+1) dimensions we 
consider a model described by the Lagrangian density 

L = L kin + L S k + L™° ass , (1) 

where L^in and L$k are the standard kinetic and Skyrme terms for the effective chiral field 
theory |],|2| and so are given by: 

L k in = ^Tr(UU + ), (2) 
L Sk = ^Tr\dJJU + , d v UU + } 2 . (3) 
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Here U is an 577(2) 2 x 2 matrix, U — $o + 2T$ and $o an d ^ are rea l scalar fields, which 
satisfy the constraint 

$2 + $ 2 = i. (4) 

The third term, the mass term, violates chiral symmetry and is very nonunique, as it is 
determined only in the limit of small pion fields. In our case we use the following form of 
the mass term: 

LZts = ^mlF^Tr[(U + U + -2)(U + U + + 2)]. (5) 

Our choice (|5]) of the mass term gives us two different absolute minima for the energy of 
the field configuration at constant U = ±1 (vacua states of the theory). Note that the more 
usual mass term leads to the existence of only one vacuum state of the theory (at U = 1): 

L m ass = \mlFlTr{U + U + - 2). (6) 

Note also that when the field configuration differs slightly from the vacuum U = 1, i.e. 
\U — 1| <C 1, the two mass terms (|5]) and (||) are very similar. 

Both the modified (|J) and the standard (^) versions of the Lagrangian possess solutions 
of the Skyrmion type. However, the modified version of the model possesses also solutions 
in the form of domain walls. The explicit expression for this solution can be easily obtained 
after a further change of variables from the ($, $o) field variables. Namely, we replace the 
($,$o) by (f,g,h) defined by 

(<f>, $ ) — (sin / sin 9 sin h, sin / sin g cos h, sin fcosg, cos/), (7) 

where / = f(r, t), g = g(f, t) and h = h(f, t),fe [0, tt], g e [0,ir],he [0, 2tt] . Note that the 
constraint (^) is automatically satisfied. The Lagrangian (TO) , when written in terms of /, g 
and h is given by: 

£ = %ifi + sin2 / flj + sin2 / sin2 9 h 2 , + k sin 2 / [{f,g,f - / 2 
+ sin 2 g {{f,h,f - fl hi) + sin 2 / sin 2 g {{h,g,f - g* hi)} - 2 m 2 V(f, g)}, (8) 

where k = 2/F 2 e 2 . 
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Consider now the case when fields g and h are constant. In this case we get from fl5|) the 
following equation of motion for the field /: 




Among solutions of this equation (|9[) there is the static soliton solution 



fsoi{x) = 2arctanexp (±m n (x — xq)). 



(10) 



In the (3+1) dimensional space-time this solution looks like a domain wall which links 
the two vacuum states U = ±1 of the theory. 

In terms of the angular variable / the Skyrmion solution becomes: 



where fsk is a Skyrmion profile function satisfying the boundary conditions /sfc(0) = ir, 
/sfc(oo) = and n = r/r. The Skyrmion solutions are characterized by an integer valued 
degree of the mapping of S 3 (compactified R 3 physical space) into S 3 - isospace. The 
analytical formula for this degree of mapping [it the winding number), when written in 
terms of the fields /, g and h, is 



As it was shown in our previous paper |T(J, a (2+1) - dimensional Skyrme-like field theory 
possesses solutions which correspond to topological waves. They behave like massless wave 
packets which carry the topological charge. It would be interesting to see whether analogous 
solutions exist in a (3+l)-field theory. 

To answer this question let us seek solutions of the equation of motion for the Lagrangian 
(§) in the form of a generalized cylindrically symmetric hegehog ansatz: 



U (r ) = cos fsk (r) + irn sin f Sk (r) , 



(11) 




(12) 



where we have used Latin indices to denote spatial variables. 



III. A TOPOLOGICAL SOLITON ON THE DOMAIN WALL 



f = f(t,x,p), g = g(t,x,p), h = rup, 



(13) 
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where n is integer and we have changed Euclidean coordinates (x, y, z) to the cylindrical 
ones (x, r, </?), 

(x, y,z) — > (x, p cos <p, p sin <p) . ( 14) 

It is interesting to note that h = rup is compatible with the equations of motion (ie the 
h equation is automatically satisfied). So we have to look at the remaining equations. In 
terms of functions / and g the Lagrangian becomes: 



lift -ft' fl + sm 2 / { 9 2 t - qI - 9 2 p ) - sin 2 / sin 2 g * 
+ k sin 2 / [(/, g x - f x g t f + {ft 9 P ~ f P 9tf ~ (fx 9 P ~ f P 9,? + sin 2 g ^(f 2 ~ ft ~ f P ) (15) 
+ sin 2 / sin 2 g $ (g? - g\ - g 2 p )} - 2 m 2 V{f, g)}2rr p dp dx, (16) 

In ( p^6|) we have used an expression more general than (||) for the mass term. Our new 
mass term depends on both variables / and g. The purpose of this generalization will 
become obvious soon. The equations of motion for functions / and g, which follow from the 
Lagrangian flTB]), are: 

£p _ mr, f ™ c f r„a ^2 _ ri 2 _ -2 „^\_ Lw 2 gym 



fa - fxx -fpp-p- sin / cos / (gt -9 2 x -g 2 P - sin 2 g^) + 9f 



m. 



+ ksinf cos / [(/ ( g x - f x g t f + (f t g p - f p g t f - (f x g p - f p g x f + sin 2 g (/ 2 - / 2 - / 2 ) 



- 2 sin 2 / sin 2 g (g\ - g\ - g 2 ) ^} + k sin 2 / [(f tt (g 2 x + g 2 p ) - (f pp + (g 2 x - g 2 t ) 



- fxx(g p - g 2 ) - 2f xt g x g t - 2f pt g p g t + 2f xp g x g p + g tp (f t g P + f P gt) + gtx(ft9x + fx9t) 
- g xp {g x f P + g P fx) - gtt{f P g P + f x g x ) - g xx {ftgt - f P g P ) - g PP {ftgt - f x g x ) + ^(f x g x - ftgt) 

+ 2 sin g cos g^(f t g t - f x g x - f p g p ) + sin 2 g (f tt - f xx - f pp + 4] = 0. 



and 



~ 2 m M 9V(f,g) 



9u ~ g xx -g pp -^ + 2 cot/ (f t g t - f x g x - f p g p ) + sin g cos g ^ + — Tf dg 

+ k [ftp{gtf P + g P ft) + ftxigtfx + g x ft) - f xp {f x g P + f P g x ) - fu(g P f P + g x f x ) 
fxx(gtft — g P f P ) — f PP (gtft — gxfx) + ^f{g x f x — gtft) — {g PP + ^f){ft — ft) — gxxif 2 — ff) 

+ gttift + fp) - Zgxtfxft - 2g pt fpft + 2g xp f x f p + sin 2 / sin 2 g (g tt - g xx - g pp + 9 -f)^ (17) 
+ 4 sin / cos / sin 2 g ^(f t g t - f x g x - f P g P ) + sin 2 / sin# cos g^(g? - g 2 x - g 2 p ) 

-sinscos^(/ 2 -/ 2 -/ 2 )] = 
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It is worth emphasising here, that setting h = mp has lead to a serious simplification of 
the problem. In fact, our choice of h corresponds to using one integral of motion n, whose 
conservation reflects the internal symmetry of the solutions of the theory ([!])/. Note that by 
introducing our hegehog ansatz we have been able to reduce the problem to that of having 
to solve only two coupled equations for functions / and g in (2+1) dimensions instead of 
having to consider the coupled partial differential equations for three functions /, g and h 
in (3+1) dimensions. 

Note also that in our ansatz the topological charge becomes: 

Q = -127m J sin 2 / sin g(f p g x - f x g p )dxdp. (18) 

We note that this expression is nonzero only when n ^ 0. Thus to check whether 
the hegehog ansatz gives us a Skyrmion solution we see that we have to find solutions of 
equations (18) with the boundary conditions: 



/(-oo, p) = 0, /(+oo, p) = vr; g(x, 0) = tt, g(x, +oo) = 0. (19) 

These boundary conditions guarantee a nontrivial topological charge of the system. Note 
that there are no nontrivial solutions of fllB"!) within the class of functions with separated 
variables, i.e. of the form 

f = f(x),g = g(p). (20) 



This can be easily verified by looking at the Lagrangian (16), in which the term 

~ [fx9 P - f P 9x) 2 

clearly indicates that it is not possible to find any nontrivial solutions of the form (pOj). 

Let us look next for the stationary solutions of eqs. fllBD. Consider the expression for 
the mass term V(f, g) of the form 

VUs) = ^sin 2 /(l-5 m cos 2 s), (21) 

where < 5 < 1. Using this form of V(f,g), and performing a numerical simulation, 
we have managed to find a solution of equations (ITS) with the boundary conditions (|T9"|). 
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This solution looks like a slightly deformed 2- dimensional (in p, if variables) Skyrmion (or 
hegehog) located on the domain wall near the origin x = 0. The picture of the fields / and g 
as well as the energy density for this solution is shown in Figure 1 using the r, x coordinates. 
One sees clearly that it corresponds to a bound states of a Skyrmion and a domain wall. 
The topological charge of the configuration is 1. 

The solution is stable with respect to the 3-dimensional (r — > er) scaling transformations 
as well as with respect to 2-dimensional ones, (x, p, if) —>■ (x, Xp, if). Notice that the stability 
condition with respect to these 2-dimensional transformations is similar to the condition for 
Skyrmions in the baby-Skyrmion model in (2+1) dimensions, where a mass term is needed 
to guarantee their stability (see, e.g. [IT| and references therein). 

In the case when 5 m — > , i.e. for the case when 

V(f,g) = ^sm 2 f (22) 

any initially formed field configuration which carries a topological charge located on the 
domain wall is unstable with respect to the scale transformations in the (p, <p)-plane. This 
is why any initially formed field configuration in the case 5 m = looks like a wave packet, 
which spreads along the domain wall. This shows that in the case of the Lagrangian ([!]) 
with the mass term given by (||) nontrivial topological configurations on the domain wall 
manifest themselves as topological radial waves. Thus the Skyrmion, located outside the 
domain wall, when captured, is absorbed by the domain wall. A necessary condition for this 
capture process is an attraction between the Skyrmion and the domain wall. 



IV. ATTRACTION BETWEEN SKYRMIONS AND DOMAIN WALLS 

In this section we show that Skyrmions can be absorbed by a domain wall and form 
field configurations described in the previous section. To demonstrate this we have per- 
formed numerical simulations of the time evolution of various field configurations describing 
Skyrmions and domain walls. In our work we used the Lagrangian (|XJ) written in terms of 
three angular fields /, g and h. We were forced to consider this, a more general configuration 
than (|1"3|), as for an isolated domain wall the field h is constant but this is incompatible with 
the hegehog ansatz which describes the Skyrmion. Thus, to study their interaction we need 
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to consider the most general fields. We see that the problem of determining the interaction 
of Skyrmions and domain walls is more complicated than the problem of looking for bound 
states of Skyrmions and domain walls. 

First we have performed the numerical evolution of a Skyrmion placed at a finite distance 
from the domain wall in the first modified model (5 = 0). We have done this for various 
relative orientations of the Skyrmion with respect to the domain wall. In eah case the 
evolution was similar for every orientation that we tried: the Skyrmion was attracted by the 
wall to finally merge with it. The structure formed by this process was then dissipated by 
radiating concentric waves on the wall. Figure 2 shows the evolution of a cross section in the 
x — y plane (the plane perpandicular to the wall) of the energy density for this process. One 
sees that, after beeing absorbed, the Skyrmion is transformed into a wave which happens to 
be topologically non trivial. 

Next we have performed the same simulation when S = 0.3. Again the Skyrmion was 
absorbed by the wall except that this time instead of being radiated along the wall, it formed 
a bound state with the wall. In Figure 3 we show the evolution of a cross section in the x — y 
plane of the energy density for this process. This configuration caries one unit of topological 
charge and is is very similar to the topological structure described in the previous section. 

We then looked at the interaction between a 2-Skyrmion configuration and the wall using 
5 = 0.3 again. As for the single Skyrmion, the two Skyrmions where attracted by the wall 
and merged with it to form a static ring-like structure on the wall. In Figure 4, we show the 
cross section of the energy density in the x — y plane and in the y — z plane (on the wall) 
after the structure had settled down. 

V. SUMMARY AND FINAL COMMENTS 

We have established the existence of bound states of Skyrmions and domain walls. De- 
pending on the form of the mass term in the theory this solution is either stable or unstable. 
In the case of 0(4)-symmetry broken down to 0(3) x Z 2 this solution is unstable and its 
evolution looks like an evolution of a wave packet on the domain wall which carries the topo- 
logical charge. We have also shown that Skyrmions can be absorbed by the domain wall and 
that their field configuration is transmuted into topological waves on the wall. Considering 
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the topological charge as the baryonic charge, this picture corresponds to the absorption of 
baryons by the domain walls. Note that, because of the different vacua to the left and to 
the right of the domain wall, Skyrmions and antiskyrmions (or baryons and antibaryons) 
behave differently. We end with the speculation that these arguments, ie the separation of 
Skyrmions and antiskyrmions by the domain wall and the observation of the absorbtion of 
Skyrmions by the walls may be relevant to the discussion of the observed baryon-antibaryon 
asymmetry of the Universe. 
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FIGURES 





FIG. 1. Skyrmion on the Wall (5 = 0.3). All expressions are shown using the r, x coordinates: 
a) /; b) g; c) Energy density 
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FIG. 2. Skyrmion absorption by the Wall (5 = 0). Energy density cross section perpendicular 
to the wall at: a) t=0; b) t=37.5; c) t=50; d) t=65 
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FIG. 4. Two Skyrmions on the Wall (5 = 0.3) : cross section of the energy density a) perpen- 
dicular to the wall b) parallel to the wall 
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